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Abstract 

Mean field games were introduced by J-M.Lasry and P-L. Lions in the mathematical 
community, and independently by M. Huang and co-workers in the engineering community, 
to deal with optimization problems when the number of agents becomes very large. In this 
article we study in detail a particular example called the “seminar problem” introduced by 
O.Gueant, J-M Lasry, and P-L. Lions in 2010. This model contains the main ingredients of 
any mean field game but has the particular feature that all agent are coupled only through 
a simple random event (the seminar starting time) that they all contribute to form. In the 
mean field limit, this event becomes deterministic and its value can be fixed through a self 
consistent procedure. This allows for a rather thorough understanding of the solutions of 
the problem, through both exact results and a detailed analysis of various limiting regimes. 
For a sensible class of initial configurations, distinct behaviors can be associated to different 
domains in the parameter space . For this reason, the “seminar problem” appears to be an 
interesting toy model on which both intuition and technical approaches can be tested as a 
preliminary study toward more complex mean field game models. 


1 



I. INTRODUCTION 


Many problems in different fields deal with a situation where many identical and 
interacting agents try to minimize a cost through the choice of a strategy. One can 
think of economic agents trying to maximize their profits, of people in a crowd trying 
to minimize their discomfort or to particles in a fluid “trying” to minimize their 
energy. 

A general framework making possible to model a large class of such problems 
has been introduced in 2006 by Lasry and Lions BIZI and Huang et al. [3] under 
the general terminology of “mean held game theory”. Largely inspired by statistical 
physics, this approach addresses the limit where the agents face a continuum of 
choices (states) in which they can evolve only locally, and the number of agent is large 
enough that self averaging processes are at work. This approach leads to a system of 
partial differential equations coupling the density of players and the optimization 
part of the problem. 

Mean held game theory has been intensively studied in the past few years, and 
in spite of its relative youth, a very large number of results have been obtained in 
the mathematical |4HSJ and socio-economic community [SHIS]. A recent overview is 
given by Gomes and Saude in [IT] . Most of the focus however has been put either on 
the conditions required to prove rigorously the existence and unicity of the solutions 
of the equations of mean held game theory [8], or on the study of particular models 
based primarily on numerical treatments [7|. A more “qualitative” understanding of 
the behavior of the solutions, based on the identification of the relevant time and 
length scales, and on the analytical study of the solution in various limiting regime, 
has received significantly less attention. 

Our goal in this paper is to perform this program for a simple model, introduced 
by Gueant et al. in 2010 [15] . called the seminar problem to be described in more 
details below. The essential point here is that this “mean held game model” is in 
some sense very close to the everyday “Physicists’ mean held” since all agents are 
interacting only through a very simple “held” which is actually a simple number, the 
time T at which the seminar actually starts. This particular feature allows for an 
analytical approach, similar in spirit to the physicists’ one: For hxed T, the behavior 
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of each agent becomes independent on the other, making the associated problem 
to be solvable to a large extent; then, for a given distribution of agents, the actual 
value of T can be evaluated by a self-consistency procedure. The main interest in 
this model is to provide a fully understandable toy model on which one can develop 
its own intuition and tools before tackling the full complexity of mean held game 
models. 

The paper is organized as follows : In Section [TT] we introduce the seminar 
problem in details and show that its resolution involves two essentially independent 
parts : a system of coupled (Hamilton-Jacobi-Bellman and Kolmogorov) differential 
equations on one hand, and a self-consistency problem on the other. Sections |III| and 


IV address the Hamilton-Jacobi-Bellman and Kolmogorov equations, respectively. 


Various limiting regimes are studied in details for both. Moreover, we show that 
an exact solutions to these coupled differential equations can actually be given in a 
closed form. The self-consistency condition determining the effective beginning of the 
seminar T is discussed in Section [VJ eventually leading to the construction of a “phase 
diagram” for this toy model. Concluding remarks are gathered in Section |VT| The 
paper is completed by three Appendices where technical computations are shown. 


II. THE SEMINAR PROBLEM 

The model 

Consider a corridor at the end of which is a seminar room. A seminar is planned 
at time t but people know that in practice, it will only begin when a large enough 
proportion of the lab members 6 (known), will be seated. 

The members of the laboratory thus move according to the following considerations: 
They do not want, to arrive too early in the seminar room because they do not 
particularly enjoy waiting idly as the room fills. On the other hand they are aware 
that the lab director and the seminar organizers will already be in the room at time 
t, and will frown upon late comers. Furthermore they really want to understand the 
content of the seminar and are concerned that missing the actual beginning might 
not help in this respect. 

For every agent, this is summarized by the following cost function associated with 
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the arrival time t : 


c(t ) = a[t-t} + + P[t - T] + + 7 [T - t]+ , (1) 

where T is the effective beginning time of the seminar. In Eq. (|T|, a, (5 and 7 are 
positive real numbers and respectively quantify the sensitivity to social pressure, the 
desire not to miss the beginning of the seminar, and the reluctance to useless waiting. 
We assume these parameters to be the same for all members of the laboratory. We 
also assume (7 < a) so that the cost c(t) is actually minimal for the official starting 
time t. 

The corridor is represented by the negative half-line line ] — 00 , 0], and the seminar 
room is located at x — 0 . At time t — 0, people leave their office to go to the 
seminar. Each member of the laboratory i — 1... N, controls her drift aj(f) toward 
the seminar room but is subject to random perturbation (stopping to discuss with 
somebody, going back to take a pen and then giving up the idea, or speeding up to 
catch up a friend for example), modelled by a Gaussian white noise of variance cr 2 . 
A given participant thus moves according to a noisy dynamics: 

d Xi = a,i(t) d t + adWi(t) (2) 

where, 


Xi(t) is the agent position at time t , 
a,i(t) is her drift at the same time , 
dllqft) is a normal white noise . 


Again, except for their initial positions, all agents have the same characteristics. 

In addition to the cost c(t) associated to the arrival time (Eq. 0 ), agents dislike 
having to rush on their way to the seminar room and the total cost function therefore 
includes a terms quadratic in the (controlled) drift a;(t). An agent leaving her 
workplace Xo at t = 0 has thus to adapt her drift in order to minimize the expected 
cost 


Jr\a] — E 


c r + 


2 Jo 


a ?( r ) dr 


associated with Eq. (|2]) and the initial condition X(t = 0) = Xo<0. In Eq. 
the first passage time at x = 0 


( 3 ) 

f is 


r = inf: X(t) = 0 } , 


( 4 ) 
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and E is the expectation with respect to the noise. 

We define N(t) as the cumulative distribution of arriving time (the percentage 
of people arrived before t). If the quorum is met before the official time t of the 
seminar, this latter starts exaclty at t. If on the other hand the quorum is met at a 
later time, T is determined by the self-consistency relation N(T ) = 6 (more formally 
T = inf{t > t : N(t) > 0}).|l6j 

Within the mean field approximation, the total number of researchers in the lab is 
assumed to be large enough that the individual choices of a given agent, and thus her 
arrival time, cannot have any significant impact on T. Each agent should thus solve 
the optimization problem Eqs. (|2])-([3]) for herself, assuming T fixed. Introducing the 
value function 


u{x,t) = min |e c(f) + ^J t a-(r) dr j 


(5) 


subject to the initial condition X(t) = x, this optimization problem is equivalent to 
the Hamilton-Jacobi-Bellman (HJB) equation (see e.g. [IT]): 


du l f du\ a 2 d 2 u 
dt 2\dx) 2 dx 2 


u(x — 0, t) — c(t ) 


( 6 ) 


and the optimal drift is given by 


a(x, t ) = —d x u(x, t ) . 


(7) 


The second hypothesis underlying the mean field approximation is that, beyond 
the total size of the agent population, the agent density itself is large enough that 
the distribution of agents is self-averaging: therefore everything happens as if at any 
given location and time, each realization of the noise was experienced by somebody. 
Assuming a (normalized) initial density of participants mo(.) at time t — 0, Eq. ([2]) 
thus implies that this density will evolve under the Kolmogorov equation (see e.g. 

my- 

dm dam a 2 d 2 m 
dt + dx 2 dx 2 

' m{x — 0, t) — 0 • (®) 

m(x,t = 0) = m 0 (x) 
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Once this equation is solved, the quorum condition 


N(T ) 



r ° 

/ m(x,T) 

J — OO 


0 


(if T > t) 


(9) 


(or N(T) > 6 if T = t) provides a self-consistent condition which has to be fulfilled if 
T is indeed the actual starting time of the seminar. 


General strategy 

The toy model that we just described depends on a few parameters which play 
different roles. Some of these parameter are just “numbers”. For instance the official 
time of the seminar t, which mainly fixes a time scale. Or the parameters a, /3 and 7 
of the cost function c(t) Eq. (jlj which, as we shall see, govern the typical amplitude 
of the drift velocity. In the same way the noise strength a will govern the diffusion 
velocity. 

Another parameter of the problem is the initial distribution of agents mo(x). Being 
a function rather than just a number it is a little bit more difficult to characterize 
simply. It defines a mean initial position (x)o, but also moments of arbitrary order, 
which may introduce various length scales into the problem. 

We are helped here by the linear character of the Kolmogorov equation. Indeed 
introducing the elementary solutions G(x,t\xo) which are the solutions of Eq. ([8]) 
with a Dirac mass 6(x — Xo) as initial condition, the solution for an arbitrary mo(x) 
is obtained through the convolution 

[0 

m{x,t) — / dx 0 G(x, t\x 0 )m 0 (x 0 ) . 

J — OO 

Therefore, introducing 

p(x 0 ,t)= f dxG(x,t\x 0 ) , (10) 

J — OO 

which thus measure the proportion of agents starting from xq who have not yet 
reached the seminar room at time t, the self-consistent condition Eq. @ reads 

[ dx 0 p(x 0 ,T)m 0 (x 0 ) = 6 (if T > t) , (11) 

J — OO 

with 6 = (1 — 6) the proportion of agents still in the corridor when the quorum is 
met. 
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The solution of the self consistent problem can therefore be split quite neatly 
in two distinct parts. The first part will be to analyze, and solve, the Hamilton- 
Jacobi-Bellman and Kolmogorov equations ([6| and ([8]) assuming T known. More 
specifically, the goal in this first part will be to compute the function p(x o, T) for 
arbitrary xq and T. This is what we shall do in the two following sections. For this 
part we obviously do not need to specify what is mo(xo). 

Once p{x o, T) is known, the self consistent problem reduces to Eq. 0 - It then of 
course involves the initial density m 0 (x o), as well as p(x o, T), but this latter quantity 
summaries all the required information, beyond mo(io)) about the system. This 
second aspect of the problem will be addressed in section |V| 


III. RESOLUTION OF THE HAMILTON-JACOBI-BELLMAN EQUATION 

Except for its rather non-standard boundary conditions, the HJB equation Eq. 0 
is closely related to a Burger’s equation, and it can in the same way be solved exactly 
through a rather standard Cole-Hopf transformation. Before we do so however, we 
find it useful to consider first the limiting behaviors of very small and very large cr’s. 


A. Small a 

To understand the regime of very weak noise, let us consider the noiseless limit of 
Eq. ([6]), which takes the form of the Hamilton-Jacobi equation 

L(d t u, d x u) = 0 , (12) 

associated with the free propagation Hamiltonian 

2 

L(E, P )=E- P - (13) 

complemented with the boundary conditions 

u(x — 0,t) — c(t). (14) 

Introducing a fictitious time £ and noting ( ) = d( )/d£ the corresponding time 
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derivative, the Hamilton dynamics associated with L is given by USED! 

i = fl = 1 E = - d -± = 0 
dE dt 

dL c)L 

x=- = -p p=-— = 0 . 
op ox 


( 15 ) 


Solution of the Hamilton-Jacobi equation are typically obtained through me 
method of characteristics. Here, this amounts to build a one parameter family of rays 
iy(£) (r = (E,t,p,x)) indexed by f, such that r f (£) is solution of the Hamilton’s 


equations Eq. (15), and with initial conditions r(£ = 0) = (E 0 ,t 0 ,p 0 ,Xo) imposed by 


Eq. (14) as 


tn — T 


Xo = o 


„ dc, . . 

S„=-(r ) ^c(r ) 


L(E 0 ,p 0 ) = 0 . 


This last equation imposes 

Po = . (16) 

Once this one parameter family of ray is build, the solution of the Hamilton-Jacobi 
equation just reads (for t < f, and thus negative £) 

u(x=x f (£),t=tf(£)) = c(f) + f (Et+px)d£ 

Jo 

= c(f) — (f — t)c'(f) — x\J2c'(t) . (17) 


As illustrated on 
different regions, 


Fig.g the quarter plan {x < 0 ,£ > 0 ) has to be divided in four 


Region (0) : 
Region (1) : 
Region (2) : 
Region (3) : 


x ^ -yj2(a + (3)(T - t) 

~\j 2 (a + P)(T - t) ^ x ^ -^/ 2 (a- 7 )(T - t) 
\j2(a--i)(T - t) < x < -^/ 2 (a - 7 )(t-t) 

— \j2 (a — y)(t — t) ^ x < 0 , 


(18) 


for which the application of Eq. (17) is somewhat different. In region (0) for instance 
the relevant rays reach x — 0 at f > T which corresponds to c'(f) = (a + f3) = c' 0 . 
In the same way for region (2) t < t < T and c'(f) = (a — 7 ) = c' 2 . 
















Figure 1. Regions of the space (x,t) (x < 0,t > 0) and the corresponding drift fields for 


o i 0 (cf Eq. (19)) [a 0 = \/2(a + /3),a 2 = sj2{a - 7 )]. 


Region (1) then corresponds to f = T, where c'{ f) is discontinuous. It can be 
easily justified (e.g. by viewing c(t) as the limit of a family of differentiable functions) 
that the correct procedure here is to use all the rays emerging from (.x — 0,t — T ) 
with all possible values of c'(f) within the interval ](a — 7 ), (a + (3)[ (and thus all 
velocities p 0 within ] \j2(a — 7 ), \j2(a + /?)[). 

In the same way region (2) corresponds to f = t, and one should use all the rays 
emerging from (x = 0,t = t) with all possible values of c'(r) within ] 0 , (a — q)[ (and 
thus all velocities p 0 within ]0, \J 2 (a — y)[). 


Note that because c'(r) is negative for t < t, Eq. (16) has not real solution for p 0 


and it is not possible to fulfil the boundary conditions Ecp (14) in this time interval 


for the Hamilton-Jacobi equation. In a small layer near the line, (x = 0, 0 < t < t), 
the fact the Hamilton-Jacobi equation is first order when the HJB equation is second 
order implies a qualitative difference between the limit of small cr’s and <7 = 0. 
Acknowledging this, and keeping in mind the procedure explained above to handle 


the discontinuities of c'(r), Eq. (17) gives an explicit solution for the Hamilton-Jacobi 
equation Eq. (12)- What will be needed though as an input for the Kolmogorov 
equation Eq. (j 8 j) is not so much u(x,t ) itself than its spatial derivative — d x u which 
through Eq. ([7]) specifies the drift a(x,t) in Eq. Q. From Eq. (17) we see that — d x u 
is just the velocity — po of the free motion on the corresponding ray. We obtain 
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therefore the following results, 


y^2 (a + f3) = ao 


in Region (0) 
in Region (1) 


—x 



in Region (2) 
in Region (3) 


(19) 


—X 


At-t) 


valid as a —)■ 0. On Fig. [lj this velocity is shown as the inverse slope of the arrows. 

In the deterministic limit considered in this subsection a finite fraction of the 
agents (namely all those starting in region (1)) arrive exactly at time T. The quorum 
condition Eq. ([£]) is therefore ill-defined in this limit. 

B. Large a 

Let us consider now the HJB equation Eq. (|6]) in the limit of very large a' s. This 
amounts here to neglect the nonlinear term |( d x u ) 2 , yielding the backward diffusion 


equation 



( 20 ) 


In this subsection it will be convenient to use for the boundary conditions a slightly 
modified version c\(t ) of the cost function Eq. (JTj) , 


c A (t) = c(t) for t < A 


( 21 ) 


c A (t) = c(A) for t > A 


with (A t, T) a very large time (one may imagine for instance that once the 
seminar is over, there is less marginal incentive to reach the seminar room). 

There are many ways to derive a solution of Eq. ([6]) , but a relatively transparent 
one consists in going back to the original optimization problem, i.e. to define u(x,t) 
as Eq. ([5]). Indeed, in the limit of very large cr’s, this optimization is straightforward: 
if the motion is overwhelmingly dominated by the noise, the best strategy for an 
agent is just to renounce paying the cost of the drift, and hope that the diffusive 
motion will bring her in time for the seminar. 
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Let us note 


1 


G 0 (x,t ) = 


: exp 


x 


( 22 ) 


V&T'l 2aH, 

the elementary solution of the free diffusion problem. The distribution of time of 
first passage at x = 0 for free diffusion started at t 0 in x 0 < 0 is given by mi 

d r ~ x o 

P(t) = — 7 L / dxG 0 (x,t-t 0 ) 

CUC J -(“tco 

The value function u(x,t ) is just the average of the cost function c(t) = C\(t) with 
this first passage time distribution. It therefore reads 

roo 

lim u(xn,tn) = / dtc(t)P(t) 

< 7->00 J tQ 

1 , c(t + to) 

= -x 0 / dt -- 

J 0 t 

[The fact that Go(x,t) is the elementary solution of the diffusion equation has been 


-Go(xo, t ) , 


(23) 

(24) 


used to transform (23) into (24) 


Note that Eq. (24) would be valid for any choice of the final cost function c(t) as 


long as the integral converges in +oo, i.e. as long as c(t) grows less than linearly at 
infinity. Thus the need to modify the large t behavior of c(t) in this subsection. 

C. Arbitrary a 


As was mentioned at the beginning of this section, the HJB equation can actually 
be solved for arbitrary values of cr. Indeed, using the Cole-Hopf transformation i.e. 
setting u(x,t) = — a 2 In cp(x, t) yields a linear equation for (ft(x,t): 

<90 | cr 2 <9 2 0 

(25) 


dt 2 dx 2 

e(t) 

(j)(x = 0 ,t) = e~~^~ 


It’s solution is thus the same as Eq. (24) with c{t) = exp [— c(t)/cr 2 } as cost function. 
We obtain in this way 

c(t+r) 


0(x, t) — — X 


T 


-Gq(x, t) dr 


u(x, t) = — a 2 In <p(x, t ) . 


(26) 

(27) 


An explicit expression of 0(x, t) in terms of elementary functions is given in Ap- 
pendix[A| (see Eq. (A4)). We just stress here that, because for t larger than T the 
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cost function Eq. ([I]) becomes linear (c(t) = (a + (3)t — (at + /3T)), 4>(x,t) takes a 


particularly simple form (see Eq. (A5)) from which the value function is deduced as 


u(x,t>T ) = —\J2(a + /3)x — c(t ) . 


(28) 


As a consequence, for times beyond T the drift a(x,t) is just the constant 


a(x,t>T ) = \j2(a + j3 ) = ao 


(29) 


From Eqs. (26)-(27), the limiting behaviors Eqs. (19)-(24) can be recovered. This 


is particularly simple, for instance, in the large cr regime if one uses the regularized 
version c\(t ) of the cost function. Indeed in that cas we see that as soon as a 2 c(A), 


we can expand the exponential in Eqs. (26) and the logarithm in Eqs. (27), and 


u(x,t) reduces to (24) in lowest order in 1 /cr 2 . Things are slightly trickier for the true 


(non regularized) cost function c(t) since however large cr maybe, the +oo limit of 
the integral in Eqs. (26) is such that c(t + t 0 ) cr 2 (which actually simply provides 


an effective cutoff for the integral). We find in this case (see appendix [B]) that 

|x|y / 2 (a+^)-c(t) 


4>(x, t ) = exp 


a* 


+ 0(a~ 2 ), 


(30) 


implying a(x, t ) = \j2(a + f3). (If one consider furthermore the diffusive regime 
do|x| cr 2 , where the above approximation is most useful, one can further show that 
Eq. (|30j) is valid up to 0(cr -3 ) corrections.) 


When a 2 —> 0, the integral in Eq. (26) can be approximated using the steepest 


descent approximation in regions (0) and (2); or noting that it is dominated by the 
boundary contributions at t + t 0 = T in regions (2) or t + to — t in region (3) (see 
Eq. (18) or Fig. [I]). Details of the calculations and the precise condition under which 
the approximation applies are given in appendix [Cj The drif velocity which can be 
expressed as a(x,t ) = a 2 d x cj)/(l) is computed along the same lines, and one recovers in 


this way exactly Eq. (19), except for a small region near {x = 0; 0 ^ t ^ t} scaling 


as a 2 in region (0) and (2) and as a in region (1) and (3). 

To conclude this section, we note that the drift a(x,t ) obtained from the agents’ 
optimisation is exactly \j2(a + (3) for t > T. but is actually close to this value in the 
entire region (0) already for small cr’s. As a’ increases, the part of the domain for 
which a(x,t) — \j2(a + j3) increases beyond region (0), and extends to essentially all 


12 






















positions and times for large cr’s. This evolution of the drift with cr is illustrated in 

Fig. m 

IV. RESOLUTION OF THE KOLMOGOROV EQUATION 

We turn now to the resolution of the Kolmogorov equation. As we shall see 
below, this equation can, in the particular case we consider here, also be solved 
exactly. Before proceeding to the description of this exact solution, we find it useful 
nevertheless once again to discuss briefly the two limiting cases. 


A. Limiting cases 

1. Case cr 2 —>■ 0 


If we just set a to 0, the Kolgomorov equation reduces to: 

dm d(a(x,t)m ) 

dt dx 

m(x, t = 0) = m 0 (x) 


= 0 


(31) 


with the velocity a(x,t ) given by Eq. (19). 

Noting the total derivative attached to the flow a(x, t), Eq. (31) reads Dm/Dt = 
0 in regions (0) and (2) of Fig. [TJ Dm/Dt = —x/{T — t) in region (1), and Dm/Dt = 
—x/{t — t) in region (3). This yields (see e.g. [1$]) 


m(x, t) = < 


m 0 (x - 2(a + f3)t} 

in Region (0) 

Ut™o (x£i) 

in Region (1) 

u 

1 

1 

o 

in Region (2) 

th m o (; x ih) 

in Region (3) 


(32) 


In other words, all agents starting at t — 0 from a position xq < — y 2(a + j3)T will 
arrive after T, all agents starting at t — 0 from a position Xq > —^2(a — y)T will 
arrive before T, and all agents between — ^j2(a~+~^)T and — ^2{a-^y)T will arrive 
exactly at time T. Therefore, at a = 0, the function p(x 0 ,t) needed to define the self 
consistent condition Eq. GD become singular at t = T. We shall see below how this 
behavior is regularized for a small but finite a. 
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Figure 2. Evolution of the drift with a in heat representation. Toprcr = 0.5, middle: a = 2, 
bottom: a = 7 [a = 2, j3 = 1, 7 = 1, t = 1, T = 2] 
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2. Case a 2 —»• +oo 


When a 2 —> +oo, we have seen that the drift velocity a(x,t ) tends toward the 

constant value a o = \j2(a + /3). Intuitively, this is due to the fact that when 

optimizing her drift, the agent does so not so much having in minds the median 

arriving time but rather try to ensure herself against possible late arrival due to the 

noise. In this limit the motion of the agents become extremely simple, as the mass of 

participants is transported by an advection-diffusion equation with constant drift a 0 

2 

and diffusion coefficient Forgetting for now the small technicalities existing for 
small |x|, this would imply that we should simply consider the Kolmogorov equation 
with constant drift 

dm dm a 2 d 2 m 
dt a ° dx 2 dx 2 

' m(x — 0, t) — 0 • (33) 

m(x, t = 0) = mo(x) 


Explicit solutions of Eq. (33) are well known ([ 21 ]), and in particular the elementary 


solution for an initial distribution mo(x) — S(x — xq) is given by 


G CD (x,t\x 0 ) = 


a 


\J‘lnt \ 


< exp 


(x — Xq — aot)' 

2 cr 2 t 


exp 


/2a 0 x\ 

V u 2 ) 


exp 


(X + Xq + a 0 t ) 2 

2 a 2 t 


(34) 

We shall see below that this expression indeed provides the leading large a asymptotic 
approximation of the true solution. 


B. Full resolution of the coupled problem 


We turn now to the solution of the Kolmogorov equation for an arbitrary a. For 
this purpose, let us write the agent density as 


m 


(x,t) = (x,t) , 


(35) 


with u(x,t ) the solution of the (HJB) equation Eq. ([6]), which is thus such that 
— d x u = a(x,t). Inserting Eq. ( [35] ) into Eq. ([8]), we find that: 

2 0 u 4 2 t, _~ (du 1 f du\ 2 a 2 d 2 u 

^,r- y Cr = r --- - + -- 


(36) 
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But vl(x 5 is ct solution of th.6 (HJB) ocjucttion Eq. The right hand side of 
Eq. (36) is therefore uniformly zero, and this equation can be written as a simple 


diffusion equation without drift 



cr 


0 2 T 


2 dx 2 


= 0 


r(x = 0,t) = 0 and r(r,f = 0) = e 


u(x,t= 0 ) 

— m — <= 


m 0 (x) 


Noting 

Go bs (x, t\x 0 ) = {G 0 (x, t\x 0 ) - G 0 (x, t\ - x 0 )) (37) 

the elementary solution (Green’s function) for the diffusion without drift but with 
absorbing boundary in zero, obtained straightforwardly using the method of image 


from the elementary solution of the free diffusion problem Eq. (22), we find 
T(x,t)= [° G- hs (x,t\x 0 )e u ^^ 2 m 0 (x 0 )dx 0 . 

■ 7—00 


m 0 (x 0 ) dx 0 , 


(38) 


Inserting the expression Eq. (27) of u(x,t ) yields 

, Gn hs (x, tlxo) 

m(x,t) = <Kx,t)J_ go 0( ^ t = Q) 

and in particular, setting the initial distribution as a Dirac mass located in xq, we 
get for the elementary solution 


G(x,t\x„) = * Go”OMl*o) 


(39) 


The self consistence equation Eq. (11) only involves p(xo,t — T), and thus we need 
to compute G(x,t\xo) at t — T, which is in the range for which (j) can be expressed 
through Eq. (|A5|). After integration over the final position x we obtain 


p(xo,T) = 


d -c 0 (T)/<j 2 


dxe xao/(y2 Gf s {x,T\x 0 ) . 


(40) 


<j)(x 0 ,t = 0) J -c 

With Eq. ( |jO| we actually obtain an exact solution of the first part of the program 
defined in section |H} Indeed, both (j)(xo,t = 0) (cf. Appendix [A| and the integral 
in the r.li.s. of Eq. ( |40] ) can be written explicitly in terms of the complementary 
error function erfc and elementary functions. We thus have an exact and explicit 
expression for the quantity p(x 0: T) required to address self consistency (cf Eq. 0 
and the discussion below). Before we do so however, we will consider the large and 
small cr asymptotics of this exact result; and relate them to the expressions obtained 
in section IIV A1 
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C. Asymptotic regimes 


For large cr’s, and more specificaly when the condition (B4) is met, the approxi¬ 


mation (30) can be used in Eq.(39) and we get 


(41) 


G(x,T\x 0 ) = e-^( c(T) - cm e-^ {x °- x)ao G* hs (x,T\x 0 ) 

= e-^ {c{ 0 ) - com G CD (x,T\x 0 ) . 

Up to the factor e _ 2 °^ c ^ _Co ^ 0)) _ ^ 7 ~ 2 ^ one thus recovers the elementary 


solution of the convection-diffusion Eq. (34) so that 
p(x 0 ,T) ~ f G cd (x, T\xq) dx 

J — OO 


erfc 


X 0 + q 0 r \ _ e _2a 0 x 0 /a 2 erfc ( X 0 ~ a 0 T 


V2 a 2 T 


V2a 2 T 


(42) 


As we shall see in section [V] this equation will be mainly useful in the diffusion 
regime where, beyond the condition (B4) one may assume Xq <C a 2 T. In this case 


Eq. (42) simplifies to 


p(xo,T) = 


2Xn 


V2na 2 T 


exp 


oIT 

2cr 2 


(43) 



Figure 3. Solid line: asympotic form Eq. (|~i2]) of p(xo,T) valid for large cr’s (i.e. under the 


condition (B4|) and for T such that (a$T -C Vcr 2 T); Dashed : linear behaviour Eq. (43) 
corresponding to the diffusion regime, [a = 2, /3 = 1, T = 1.1 and u 2 = 49]. 
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Another asymptotic regime is obtained when the diffusion time t a = x^/cr 2 is 

much larger than the drift time t d (xo) (defined as the arrival time for the participant 

located at x 0 when a = 0; thus here t A = |x 0 |/a 0 in regions (0) t d = \x 0 \/a 2 in regions 

(2) t d = T in region (1) and t d = t in region (3)). For small cr’s this condition applies 

for most of the Xo axis, except a boundary layer in the region xi < x < 0, where Xi 
x 2 

is defined by cr 2 3> 

t 

In this asymptotic regime the Laplace method ( |23] ) can be used to evaluate the 
integrals occuring in Eq. (40) (we detail the calculation of 0(xo,O), in appendix [Cj 
and the evaluation of the numerator is done along the same lines). We get 


p( x o, T) = < 


1 , if xq < —Tao 

n 2 

if — Ta 0 < x 0 < —Ta 2 

«2 

0 , if x 0 > —Ta 2 


J' 2 ^2 

a 0 ~ a 2 


(44) 


We see in this way how the singular behaviour of the strict a = 0 limit is regularized 


for small but non-zero sigma (cf Eq. (32) and the discussion below). An illustration 
of this function is given on Fig. |4} 



Figure 4. Asympotic form of p(x d ,T ) for a small. Full line: T > t\ dashed: T = t. 
[a = 2, (5 = 1, 7 = 1 , T = 2], 


As discussed in appendix [Cj Eq. (44) recpiires to be valid that T is sufficiently 
distant from t (cf Eq.(C6) for the precise condition). If for instance T — t one would 
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have instead 


p(xo,T=i)={ J __ c / 3 
c'o - c ' 3 


if Xq < —ta 0 
if — too < Xq < 0 


(45) 


with c' 0 = (« + (3 ) and c ' 3 = —7 the slopes of c(r) for r > T and for r < t respectively. 
The shape of this function is illustrated as a dashed line in Fig. |dj Note that Eq. (45) 
is valid only for |&o| large enough that the motion is dominated by convection, 
and in particular does not apply at a ; 0 = 0 , where in any case one should have 


p(x = 0,t) = 0. The condition of validity of Eq. (45) can be shown for small ex’s to 
read |xo| 3> (aot 3 // 2 /a) exp(—yt/cr 2 ). 


V. SELF-CONSISTENT CONDITION 


It is time now to answer the question: “when does the meeting start ?”. Answering 
this question implies taking into account the coupling between agents mediated by the 
mean-held condition, and means in practice solving the self-consistent equation (0 


with the form of the function p(xo,T ) given by Eq. (40). 

One thing worth noticing already is that having an explicit expression for p(x 0 , T) 
could provide an alternative route - to the one given by Gueant et al. [13] - for the 
proof of the existence of the solution for T, which is associated with the continuity of 
p(xo,T). This route is of course restricted to particular models such as the present 
one. In this section however, we are not so much interested in this “proof of existence” 
than into a qualitative description, and whenever possible a more quantitative one, 
of the behavior of T as a function of the various parameters of the problem. 

A stressed in section [TTJ among these parameters the initial density of agents mo(xo) 
plays a specific role. Indeed, the other parameters, namely (a, (3, 7 ) characterizing 
the cost function c(t ), the official time of seminar t, and the intensity of the noise 
a, enter through the function p(x 0 ,T ), and their specific role has been discussed at 
length when analyzing the property of this function. The initial density m 0 (x 0 ) on 
the other hand enters only now in the discussion since the behavior of the agents is 
coupled only through T. Furthermore, mo(xo) being a function, it may have a infinite 
variety of shape, and it is clearly not realistic to discuss the more esoteric among 
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them. In the following, we shall therefore restrict our study to initial distributions 
mo(xo) that can be characterized by their mean value (xo) and their variance E 2 , 
and thus implicitly assume that E sets a scale below which the variations of mo(xo) 
are small. 


The mean value (xq) will mainly determine how much T is influenced, or not, by 
the official time of the seminar. Clearly, if (xq) is close enough to zero, almost all 
the mass will be close to the origin of the negative semiaxis and there is a point 
where the noise a will be sufficient to fill the seminar room, and the quorum will be 
met before the official beginning time, giving T — t. For larger, but not too large 
|(xo)|, agent have a real possibility to arrive in the seminar room near, or even a bit 
before, t, which will influence their optimization choices, and eventually lead to a self 
consistent T which depend on t, although T > t. For very negative (xo) on the other 
hand, there is very little chance for an agent to arrive before t. Indeed, as we have 


seen in section III, it is never optimal for a agent to choose a drift velocity higher 
than a 0 = ^2cq, where c' 0 is the slope of the cost function c(t) for time t >T. As a 
consequence, if |(xo)| 3> a 0 f, the agents determining the quorum condition (ie the 
last ones to arrive before the quorum is met) will never consider the possibility to 
arrive before t, and thus the official starting time of the seminar will play no role in 
setting T. 


The parameter E on the other hand will balance the effect of the Kolmogorov 
diffusion term in the determination of T. Indeed, a set of agents starting from an 
identical initial location will have spread on a distance ay/T at time T. So for 


E <C aVT 


(46) 


the diffusion will essentially erase any of the initial features of m 0 (x o), while for 


E S> aVT 


(47) 


diffusion plays little role for the transport of m(x : t). Keeping in mind this general 
picture, we turn now to a more detailed description of the various limiting cases. 
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A. Self consistent condition in the diffusion regime (large cr’s) 


We characterize the diffusion regime by the fact that a is large enough (condition 


(B4)) and that, for the relevant positions Xq, the time of drift is much larger than 
the diffusion time, i.e. here 

\x 0 \a 0 < cr ;. (48) 

We consider successively narrow initial distributions and wider ones. 

1. Narrow initial distributions 


A narrow initial distributions corresponds to a configuration where the initial 
width E is significantly smaller than the spreading cry/T acquired because of the noise 
during the transport - the notion of narrow initial distribution is thus e-dependent, 
and this configuration will typically be met when the noise is rather large. In that 
case the details of the initial distribution become irrelevant and mo(xo) can be 
approximated by a Dirac function 5(x 0 — (^o)))- The integral equation 0 therefore 
just become the simple standard equation 


p({xo),T) = 6 , 


(49) 


in which mo(xo ) is entirely characterized by its mean value (xq). In the large noise 


regime the conditions (B4) and (48) hold. Furthermore, we will see that for the self 


consistent value of T obtained at the end of the process one gets 


(*o> 2 « o 2 T ;. 


(50) 


Under these conditions, we can use the approximation Eq. (43) for p(xo,T), and 


Eq. (49) reads e u /y/mi = #cr 2 /(ao((ro)) 2 , with u = OqT /2 


cr- 


If 9 ao|(xo)|/u 2 , which deep in the diffusive regime will usually holds except 

for very small 9, we get in leading 1/cr order, T = sup (T*,t), with 

2 (x 0 > 2 


T* ~ 


7T a 2 6 2 


(51) 


(If T = T*, the condition (50) then just amounts to have 9 • Cl, which we assume. If 
T = t, the condition (50) is even more easily fulfilled.) T* is proportional to (xq) 2 , 
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and we recover the intuitive result that if the initial distribution is located too close 
from the seminar room, the noise fills this latter before the official starting time, 
giving T — t. 


For very small 6 , there is a range of cr 2 for which even in the diffusive regime (48) 
one has 6 <C ao|(xo)|/cr 2 . In that case 


2 a 2 fa 0 
T ~ —log 


l o 


\(*o)\ \ 

y/lT(J 2 6 J 


(52) 


and one can check that (50) holds. 


2. Wide initial distributions 


When E ay/T - which since we assume here a 2 (c(0) — Co(0)) implies 
fairly large E’s ~ the convolution with a Gaussian of width ay[T barely change 
the distribution. Everything appear then as if the Kolmogorov approximation was 
dominated by convection. 

Let us introduce Xg such that 



mo(xo)dxo = 0 , 


(53) 


which is thus the position of the participant such that a fraction 6 of the agents is 
more distant from the origin. The beginning of the seminar is entirely determined by 
the time at which the agent starting from this location and evolving in a deterministic 
way under the influence of the drift a(x,t ) (ie ignoring the effect of the noise) will 
arrive. 

In the large noise limit that we consider here, the drift is constant and equal to 
a 0 , and this just gives 


T = 


Xt ? 

do 


(The fact that Xg is necessarily of the order of or larger than E, together with (B4), 
implies T > t.) 


B. Self consistent condition in the convection regime 


In the convection regime, and more precisely under the conditions (C4)-rtC6b, the 


function p(x,T ) is well approximated by Eq. (44). We consider below how Eq. 0 
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can be solved for this form of p(x, T ) for different ranges of the initial distribution’s 
width £. 


1. Narrow initial distributions 


For very narrow initial condition, Eq. (fllj) can be as before replaced by Eq. (49) 


which, with Eq. (44), is solved as 


T = 


I (so) I 

a{6) 


(54) 


a(9) = Jd 2 2 + (al- a%)9 . 


For small 9 , Eq. (54) corresponds to a{9) ~ a 2 , ie to a (xq) near the lower border of 


region(l), for which the condition (C5) might not be fulfilled. Re-inserting Eq. (54) 


into (C5) we indeed see that Eq. (54) applies only if 


a 2 « a 0 |(i 0 >|9 2 ■ 


(55) 


For larger cr’s - or smaller 9 - we need to use for f(x, 0) in Eq. (40) the uniform 


approximation Eq. (C7) valid for x near —a 2 T. Writing T = T 0 + ST with T 0 = 


\{xf)\/a 2 , we find then in linear order 


5T 

% 


7T On 


a\ 


i9 0 ~9), 


(56) 


with 6*o = \J 8cr 2 /7r| (x 0 ) |o 2 (a|/ (al - a 2 )). 

The expressions ([54|-(|56j) are clearly independent of the official beginning time of 


the seminar t. The condition (C6), which is actually required for Eq. (44) to apply, 


indeed implies that T is sufficiently above t to become independent of this latter. 

Once | ((Co) | diminishes, and more specifically when it reaches a value close to 
a 2 t or smaller, T will on the other hand approach t. It may be interesting then to 
determine under which condition one has exactly T — t, i.e. when the quorum is met 
before the official beginning time t. 


In the convection regime, p(x,T — t) is described by the expression Ecp (45), and 
the self consistent condition to obtain T — t is that 


p({x 0 ),T = t ) < 9 . 


(57) 
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The approximation Eq. (45) is however bounded from below by its value at zero, 


7 /(a + /3 + 7 ), which is typically of order one. If, as we assume, 6 is small, Eq. (57) 


will thus not have any solution for (xo) is the convection regime. As long as the 
motion of the agents is convective, they will manage to fill the seminar room after, 
though possibly barely, the official start of the seminar. 

If |(iCo)| becomes so small that the motion at such distance is dominated by 
diffusion, then again one can eventually reach a point where the quorum is met 
before t, implying T — t. For small cr’s it can be shown that this happens when 


l(*o)l = ,A/2(a2P'7<r)(9/(l -«))exp(- T (7a 2 ). 


2. Wider initial distributions 


If the width of the initial distribution is non negligible, we need to distinguish two 
cases. For intermediate values of E, namely for E’s such that once self-consistence is 
obtained most of the initial distribution is in the range ] — aoT, — a 2 T[, we can use 
that in this range the function p(x, T) ~ ( x 2 /T 2 — a|)(ag — a|) is a simple polynomial. 
The convolution with mo(x) thus simply leads to 


dx 0 p(x,T)m 0 (x) = 


((x o y + E 2 )/T 2 - a 2 
a 0 ~ a 2 


and Eq. (54) has just to be replaced by 

T = 


+ («o - a !)^ 

The constraint that the initial distribution fits within 


\ 


(x„) 2 + s 2 


(58) 

a 0 T, —a 2 T[ implies that 


Eq.(58) applies only if E -C 0|(x o )|, i.e. for not too small 6. For smaller 6, explicit 


(but less transparent) expressions can be written down under the less restrictive 
condition E < (ao — a 2 )T ~ ((do — a 2 )/a 2 )\i(Xo)\ for specific forms of the initial 
distribution (eg Gaussian). 

If now E is large not only on the scale <j\/T but also on the scale (ao—a 2 )T, another 
approach can be used. Subtracting Eq. (53) to Eq. 0 and neglecting the variation 
of m 0 (x 0 ) near xq in the whole region ( 1 ), we can write that 2 a ' 2 dx 0 = 

0 <2q a 2 

f-a oT dxo, which implies 

o /^2 „ 2 \ 

T = ~^\ x e\ 


2 r ’ l/| \ a?, — an 


As before, this results apply only if T — t is large enough for Eq. (C 6 ) to be fulfilled 
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C. “Phase diagram” of the seminar problem 



- 21 og( 0 ) 

log( 0 ) 


Figure 5. Phase diagram of the seminar problem, in the (a 2 , |(a:o)|) plane, for narrow initial 
distributions (see text for the detailed description of the various regimes). The vertical 
dashed-dotted line correspond to the vertical cuts used in Figs. [ 6 ]an d 0 For this illustration, 
the parameters of the problems have been taken as a = 2 , = 1,7 = 1 , t = 1 , 0 = 0.2 


The results of the previous subsections can be summarized into “phase diagrams” 
such as the one shown on Fig. [5] for narrow initial distributions (equivalent phase 
diagrams can be constructed in the same way for wider initial distributions). Keeping 
in mind that, except for the transition between T — t and T 7 ^ t, there is of course 
no true phase transition here, and that the lines representing the limits between 
various regimes should be though as crossover regions (thus with a finite extension), 
we can distinguish the following “phases”: 


Region I corresponds to a motion dominated by convection, and such that the initial 
distribution is far enough from the seminar room that the initial time of the 
seminar becomes irrelevant. This region is split into two subregions. In the 


first one, I a , T = |(x 0 )|/a(^) with a(6) defined by Eq. (54). In the second one, 


lb, the fact that 0 C 1 and thus that a{9) ~ 02 makes it necessary to use the 


uniform approximation Eq. (C7) for (j>({xo ),T). In that case T = \(xq)\/ a, 2 +ST 


where ST is given by Eq. (56) 


Region II corresponds to a motion dominated by diffusion, and such again that the initial 
distribution is far enough from the seminar room that the initial time of the 
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seminar is irrelevant. Region II, too, has to be divided in two subregions. In 


II a , T = T* with T* given by Eq. (51). In II b the smallness of 6 should be 


taken into account, leading to Eq. (52). 


Region III corresponds to a motion which can be dominated either by diffusion (region 
IHa) or by convection (region Illb), but such that in any case the quorum 
is met before the official beginning time of the seminar. This region thus 
correspond to the phase T — t. 


Region III corresponds finally to a configuration such that the quorum is met slightly 
after the official beginning time of the seminar, so that T is different from, but 
close to, t. 


As an illustration, we show in Fig. [6] and [7] two vertical cuts in this phase diagram, 
in which are displayed the variations of the self-consistent time T as a function of 
(xo) for two (fixed) values of the noise a, one “small” and one “large”. For the small 
diffusion coefficient case Fig. [6] we observe, as expected from the phase diagram, a 
transition between a domain where T — t and a domain where T = |(xo)|/a(d). For 
the large a case Fig. [7J we observe, again as predicted from the phase diagram, a 
richer behavior, with the same limiting behaviors for very large and very small |(x 0 )|, 
but a larger number of intermediate regimes. 


VI. CONCLUSION 


In this article, we have considered a simple toy-model, sharing many of the 
characteristic features of generic mean field games, but with the essential simplification 
that the “mean held” actually reduces to a simple number (the actual starting time 
of the seminar T). The study of this problem can then be divided in two essentially 
independent parts : on one hand the resolution, for arbitrary T, of the system of 
partial differential equations Eqs. describing the coupling between the agents 

optimization decisions and their motion; and on the other hand the self-consistent 
condition Eq. ( [TTj ) determining the value of T. 

The first part of this program can be performed essentially completely. Indeed an 


explicit expression Eq. (40) can be obtained on a very general basis for the function 
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|<x 0 >| 


Figure 6. |(xo}| dependence of the time T for narrow initial distributions and a small value 
of the noise parameter (a = 0.6); which corresponds to the left vertical dashed-dotted line 
in the phase diagram Fig. [5| The numerical labels correspond to those of the different 


regimes in Fig. |5j Full line : numerical value obtained from the exact expression Eq. (40). 
Dashed : asymptotic expressions in the corresponding regime (see text). The parameters 
of the model are the same as in Fig. [HJ except for 9 = 0.1 which, to enhance readability, 
has been slightly decreased. 


p(x o, T ) required to discuss self-consistency. From this general result transparent 


asymptotic expressions are derived in the relevant limiting regimes (cf e.g. Eqs. (42), 


(44) and (45)). From these and the self-consistent condition Eq. (11) the qualitative 
behavior of T, and explicit expressions in many limiting regimes of interest, can be 
obtained. 

We will not re-list here the various results derived for T in section [Vj but the main 
features may be summarized as follow. An important point is that the the slopes 
Cq = (a + /3) and c ' 2 = (a — 7 ) of the cost function c(t) (cf Eq. (0) can be associated 
with drift velocities ao = y/ 2 c ' 0 and a 2 = \J 2 c 2 which fixes the scale of the drift 
velocities of the problem. Together with the characteristic length l 0 associated with 
the initial distribution of agent m 0 (x 0 ) (namely the center of mass, i.e. l 0 = |(x 0 )| for 


narrow distributions, or l 0 = \xq\ (cf Eq. (53)) for wide distributions) and the value 
of the noise parameter cr, they organize the system “phase diagram” of the problem, 
and define the relevant limiting regime. For instance if |Zq| not and \Iq\ cVt 
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Figure 7. |(a;o)| dependence of the time T for narrow initial distributions and a large value 
of the noise parameter (<r = 40); which corresponds to the right vertical dashed-dotted 
line in the phase diagram Fig. [5] The numerical labels correspond to those of the different 
regimes in Fig. [5j Full line : numerical value obtained from the exact expression Eq. (40). 
Dashed : asymptotic expressions in the corresponding regime (see text). The parameters 
of the model are the same as in Fig. [ 5 ] except for 9 = 0.01. 


the system is completely dominated by convection, and T ~ lo/a, with a equal to 
ao for large noise and to some weighted average between ao or 02 for smaller noise. 
Or if l 0 -C ci 2 t one is essentially guaranteed that the quorum will be met before the 
official starting time t of the seminar, and T = t, etc.. 

One point worth being stressed however is that in most circumstances, either the 
agents start from a location close from the seminar room (/ 0 -C do,it) and T = t , or 
they are initially far from the seminar room (Z 0 aoqt) and T becomes relatively 
quickly independent of t. The transition region for which T > t but keep some t 
dependence is actually rather restricted. This notion of a effective starting time of 
the seminar which is independent of its official starting time is clearly a bit disturbing, 
especially from the viewpoint of the seminar organizer. 

This feature can be tracked back to the fact that there exist an initial time t — 0 
at which the all agents start their optimization process and motion, and what we 
see is that this initial time plays a role which is at least as important as t in the 
determination of T. One can of course imagine that this initial time has some physical 
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signification (e.g. time at which the organizers ring a bell, etc..). One could also 
modify slightly the model to remove the reference to a uniform initial time and have 
the time r 0 at which a given agent leave her office taken as a parameter entering in 
the optimization decision. Assuming the marginal cost of not being in ones office is 
described by some parameter rj. this would amount to replacing the cost function 
Eq. (§ by 


Jt [ol — E 


c(r)-r]T 0 + - af(r) dr 
z Jto 


The problem could be analyzed along the line of what we have done in this paper, 
and would lead to a stronger dependence of T in t. 

More generally, various variations of the problem can be easily studied with 


the approach followed in this paper. In particular, Eqs. (39)-Eqs. (26) are valid 
for essentially arbitrary cost functions c(f), and could be studied for instance in 
circumstances for which the self consistent condition Eq. has more than one 
solution. The seminar problem is therefore a very versatile model, and the very 
thorough understanding of its behavior obtained in this work should help develop 
the intuition on the properties of more generic mean field game models. 
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Appendix A: Exact Solution of HJB Equation 


In this appendix we derive exact expressions for 4>(x,t) (see Eq. (26)) 



M 


r dr e -M 

\/‘2'K0 2 J 

0 T 3 / 2 


0(M) = 

where c(t) is piece-wise linear. 

c(t) = a[t - i] + + /3[t - T] + + q[T - t] + 


(Al) 


(A2) 
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Using the explicit expression for c(t) and the fact that t <T, one may write 





'yi.T—t) 

dr - 

-Af_, 

— 

X 

e 

" 2 / 

Jo 

r 3 /2 e 

0-“ V 




a(t—t)+~f(T-t) 

,[T-t] + 

dr - 

+ 

X 

e 


T-t]+ 

r 3/2 C 




a(t—t)+/3(t—T) 

poo 

dr — 

+ 

X 

e 

- 2 

l[T-t]+ 

7-3/2 e 


We define the following integral m 


I(a, b, t ) 



-ar - 

Jo r 3 / 2 


b 

T 


= e 


2 y/ab 





(A3) 


where the function erfc(z) is defined on the complex plane, liin x ^ +oc erfc(x) = 0, 
lim 3 ,_ ) ._ 0O erfc(x) = 2. Accordingly, lim^ 0 I(a,b,t) = 0 and lim t ^. +00 I(a,b,t) = 
2 e —2 y/ab f or a > o and b > 0. 

Thus the function 0(x, t ) reads 


(j){x,t) = -e 

1 _ 

+ 2 6 
1 _ 
+ 2 6 




_t (Tzf) . 7 x 2 

j (__L - 

1 a 2 ’ 2 <T 2 

. aft-f)+ 7 (r-t) , _af_ 

' 2<r 2 ’ 2cr 2 ’ ^ 

a(t-t)+p(t-T) . _ ap\x\ fl 2 

<r 2 (2e <r 2 -Z( ° 


X 


2cr 2 ’ 2cr 2 


,7-i] + )), 


*]+)) ( A4 ) 


(where a 0 = \/ 2 (Vt + 7). a 2 — y/2(a - 7 )). 

Eq. (A4) applies for arbitrary values of time and position. It takes however a 
simpler form in some time intervals: 


- For t > T the expression reduces to 


c(t) 

' rr2 


4>(x, t ) = e exp 


O' 


= e 


exp 


: y / 2(a + /3)|x| 


af ) (T—t) — (to x 
2^2 


(AS) 
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For t <t <T Eq. (A4) can be written as 


e -a(T-t)/a 2 

4>{x,t) = ---x 


exp 


exp 


2 

a\(T—t) + 2a,2x' 


erfc 


-x - a 2 (T—t ) 


2cj2 v '2rr 2 (7'—/) 

'a\(T—t) — 2a 2 x\ f— x + a 2 (T—t) 


2 a 2 


erfc 


exp 


(al(T-t ) + 2 a 0 x\ ^ (x + a 0 (T — t)\ 
V 2cj2 ) er C V y/2a 2 (T-t) ) 

_ ( a 2 0 (T-t)-2a 0 x \ ( - x + a 0 (T-t ) 

p l 2 ° 2 r { 


At t — 0 Eq. (A4) reads 


1 7 T ru rp2 


cr z 2(J 2 
at+f3T . y/2( a +ff)\x I 

+ -e cr 2 ( 2 e <z 2 


(7 


2 cr 2: 


1 cr 2 ’ 2 a 2 ’ J 


(A 6 ) 


(A7) 


Finally, we stress that for large cr’s, and more precisely under the condition 
{a,/3, , y)i -C a 2 (but irrespective of the value of T and x ) Eq.(A6 ) provides a good 


approximation of the exact 00 c, 0 ) for (t < t) (t and t can then be set to zero in this 
equation). 


Appendix B: Evaluation of 0(x,i) for large u’s (diffusion regime) 


In this appendix, we evaluate the large a asymptotic of the function 00r, t ) defined 


by Eq. (26). For this purpose, let us introduce c 0 (t) = a(t — t) + f3(t — T), the linear 
f un ction such that c(t ) = c 0 (t) for t > T, and 

c o(t + t ) 


,, , _ f°° dr ( c 0 (t + r)\ t N 

00 [x,t)=-xj — exp I-— G 0 (i,r) 

= exp + aokl]) 


(Bl) 

(B2) 


(this last expression exactly corresponds to Eq. (A5), and is obtained in the same 
way). 
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The difference between <j)(x,t) and its approximation cp 0 (x,t) can be expressed as 


(M) - 0 O (M)| = -X f 

Jo 


00 dr f co(t + r) 


0 r 


exp 


cr- 


G 0 (x, r)K(t + r) (B3) 


where 


K (t) = 1 — e 


— (c(r) - c 0 (r)) 


a 


is a positive (since c(r) > Co(r)) decreasing continuous function which is uniformly 
zero for r larger than T. We thus have | (j>(x,t) — <f>o(x,t)\ < (j) 0 (x, t)K(t). 

As soon as 

a 2 > (c(t) - co{t)) , (B4) 

K(t) is 0(cr~ 2 ), and therefore 


<p(x,t) = tj>o(x,t)(l +0(cr ! )) . 


(B5) 


Eq. (B5) applies whenever the condition (B4) is fulfilled. In practice however, it 
is mainly useful if the diffusive regime, i.e. when 


ao|;c| -C u 2 


(B6) 


(ao = y/2(a + f3)). In that case, since Eq. (B4) morally implies a 2 |co(r)| (this is 
clear as soon as t < (at + f3T)/(a + f3) since then c 0 (t) < 0, but remains generally 


true unless r ~ T), one has cj) 0 (x,t) = (1 + 0(a~ 2 )), and Eq. (B5) provides little 
information on the variations of (f>(x,t). 

It may be therefore interesting in this case to compute the 0(c r~ 2 ) corrections. 
Noting that K (r) = 0 for r > T, we have 


0OM) - 0o(x,t) = -x [ 

Jo 


T-t (J r 


r 


exp 


c 0 (t + r)' 


cr z 


G 0 (x,r)(K(t)+5K(r)), (B7) 


where 5K(r) = (K(t + r) — K(t )). The term involving 5K , which is linear in r near 
0 (and thus do not benefit from the r” 1 / 2 divergence) can be shown to be 0(a -3 ) 
relative to (j> 0 (x,t), and we get 


0(M) = (0o(^,^) + K(t))( 1 + 0(c r 3 )) 

= exp ( _ ^2[ c ^) + a oMl) i 1 + , 


valid therefore when both conditions (B4) and (B6) apply. 


(B8) 

(B9) 
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Appendix C: Evaluation of 4>(x,t) for small ex’s (convection regime) 


In this appendix we evaluate the small a asymptotics of (p(x, t ) using the saddle 
point approximation (in regions (0) and (2)), and more generally the Laplace method. 
Introducing 

$(r) = c(t + r) + , 

the integral we want to compute is of the form 

$(r)\ 

a 2 J ’ 

with f(r) = — x/r 3 / 2 , and will therefore be dominated for small u’s by the minima 
of <L(t). 

The condition $ r (r*) = 0 leads to the equation 

t* = —x/\j2d(t + r*) (Cl) 

which admits a solution in regions (0) and (2) but not in regions (1) and (3). One 
therefore has to use the saddle point approximation in regions (0) and (2), and 
boundary contributions is region (1) and (3). 


(p(x,t) = 


yphx~o^ Jo 


drf{r) exp 


1. Regions (0) and (2) 


In regions (0) and (2), the stationary point is 


t* = —x/ao [in region (0)] 

t* = —x/a 2 [in region (2)] , 


which, with ^"(r) = x 2 /t 3 , gives within the saddle point approximation 


c f>(x, t ) ~ exp 
= exp 


(-v\ 

^ a(T — ^ p + 2ao,2^^ 


(C2) 


This approximation is valid as long as a 2 <C |a;|ao, 2 , or in other words as long as the 
ratio between the drift time td = |x|/ao ,2 and the diffusion time t a = x 2 /a 2 is small. 
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2. Regions (1) 


In region (1), there are no solution to Eq. (Cl) as the minima of <3>(t) correspond 
to a discontinuity of the cost functions c(t + r) (at t + r — T). Linearizing <t>(r) on 
both side of this discontinuity and neglecting the variation of /(r) we get 


, 2{T—t) 

<p{x, t) — —xcm / ---exp 


71 


X 


G 

1 


c(T) + 


x 


2 (T-t) 
1 


(C3) 


_al(T-ty-x 2 cl 2 2 (T -t) 2 - x 2 _ 

Noting a(x,t ) = x/(T — t ) the drift velocity of an agent within region (1) in the 
a —» 0 limit ( a(x,t ) G [a 2 ,a 0 ]), t d = {T — t ) the drift time and t a = x 2 /o 2 the 
diffusion times, Eq. (|C3|) applies under the condition that: 


U 

ta 

td 

ta 


a 2 (T — t) 


x “ 


a 2 (T-t) 


x* 


< 1 


< 1 - 


a 2 [x , t) 


a 2 < 


d 2 (x,t ) 
(T - 1 ) 


(d 2 (x,t) 


(C4) 

(C5) 

(C6) 


The two hrst conditions express that if generally speaking Eq. (C3) requires that 


the time of drift is much shorter than the diffusion time, the requirement becomes 
more and more stringent as (x,t) get closer from the boundaries of region (1) where 


a(x,t ) —> a 0 or a(x,t) —>■ a 2 . The last condition signal that Eq. (C3) is valid only if 
T differs significantly from t. 

The calculation of 4>(x, t ) in region (3) in the small a limit proceeds essentially 
along the same lines. 


3. Uniform approximations 


The conditions (C4)-(C5) express that the transition between regions (0) and 


(1) (i.e. x ~ — a 0 (T — t )) as well as the transition between regions (1) and (2) (i.e. 
x ~ — a 2 {T—t)) need to be treated a bit more carefully and require the use of uniform 
approximations. 

For (a, /3, -C a 2 , one way to derive these uniform approximation is just to select 


the dominating contribution of Eq. (A6). Indeed, deep in region (0) (respectively 
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deep in region(2)), one can check that the third term (respectively the first term) of 
Eq. ( |A6[ ) in which erfc is replaced by its asymptotic value 2 recovers exacly Eq. ( |C2| ). 
Near x = —ao(T—t) or x = —a 2 (T—t) the uniform approximation amounts to keep 
the full dependence of the erfc. 

For instance near x = —CL 2 T 




3 -a(T-i)/o - 2 


exp 


' a\(T—t) + 2 a 2 x' 

2 ^ 


erfc 


—x — a 2 (T—t ) 

v^cr-t) 


(C7) 


This expression will interpolate smoothly between Eq. ( |C2[ ) and Eq. ( |C3| ). [This latter 
can be seen as being obtained using the large x asymptotic erfc(x) ~ exp(— x 2 )/Vttx 2 


for the four terms of Eq. (A6) 


4. drift velocity 


With the knowledge of 4>(x,t), the drift velocities in the small sigma regime can 
be obtained from the spatial derivative of u(x,t ) = — a 2 log (j)(x, t) (cf Eqs. ([7]) and 
(j27[)). I 11 leading a order Eqs. (C2) and (C3) yield 


a(T — t) — S(T—t) — aox in Region (0) 


u(x,t ) ~ < 


c(T) + 


x 


2(T - t) 


in Region (1) 


a(T — t) — -^-(T—t) — a 2 x in Region (2) 


(C8) 


(the expression of u(x,t) for region (3) can be obtained in the same way). Taking 


the spatial derivative of these expressions yields the velocities Eq. (19). 


Alternatively, one can compute the drift velocity a(x,t) = a 2 d x (p(x,t)/(p(x,t) from 
the spatial derivative d x (f)(x,t), which can be evaluated following exactly the same 
steps as for (ft(x,t). This of course gives the same result. 
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